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Abstract. We comment on the article by Fujimoto (1997 J. Phys. A: Math. 
Gen. 30 3779) where the exact equilibrium crystal shape (ECS) in the critical 
Q-state Potts model on the square lattice was calculated, and its equivalence with 
ECS in the Ising model was established. We confirm these results, giving their 
alternative derivation applying the transformation properties of the one-particle 
dispersion relation in the six-vertex model. It is shown, that this dispersion 
relation is identical with that in the Ising model on the square lattice. 

In paper H Fujimoto determined the equilibrium crystal shape (ECS) in the Q-state 
Potts model on the square lattice at the first order transition temperature. Fujimoto 
claimed [Q, that ECS is universal for a wide class of models, including the eight- 
vertex model ^, D , and the Ising model on the square lattice Q] . The origin of this 
universality is still not well understood. 

The subject of the this comment is to show, that the one-particle spectrum in the 
six- vertex model (which is known to be equivalent to the critical Q-state Potts model), 
and in the Ising model on the square lattice are also the same. So, one can say, that 
the universality of ECS in different square-lattice models reflects the universality in 
the one-particle dispersion relation. 

Consider the Q-state Potts model (Q > 4) on the square lattice, shown in figure 
|l|. The sites of the Potts lattice are depicted by open circles. The model Hamiltonian 
is given by 

3, I J, I 

Here we use the same notations as in H, constants K\, K 2 obey the critical 
temperature condition ||: (exp K\ — l)(exp K 2 — 1) = Q. Let t/j be a vector associated 
with a horizontal row of the Potts lattice, ip = il>{.~, cr±, a 2 , ...,07, ...). Denote by T\ 
and T2 the two shift operators El |) , shown schematically in figure [l]. Operator T2 is 
the row-to-row transfer matrix of the Potts model, and T\ acts on the vector ip as 

T 1 tl)(...,a 1 ,a 2 , ...,a h ...) = ip(...,a 2 ,a 3 , ...,a t+1 , ...). (2) 

It is well known, that the critical Potts model is equivalent to the six- vertex model 
j| ^ , and to the interaction round a face (IRF) model 0, [| . Consider for defmiteness 
the associated six- vertex model. Its sites are shown by full circles in figure |l|. One can 
define two shift operators T x and T y corresponding to translations along the directions 
x, and y. Note, that operator T y is the row-to-row transfer matrix of the six-vertex 
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Figure 1. The Potts model and the associated six- vertex model square lattices. 
The Potts model sites and two-spin interactions are shown by open circles, and 
by broken lines. Indices I and j enumerate rows and columns of the Potts lattice. 
The sites of the six-vertex model are shown by black points. 



model. As it is clear from figure || these shift operators are related with T\ , T2 by the 
equations 

T\ = T x T y , T 2 =T- x T y . (3) 

Consider operator B(9), generating a one particle excitation from the ground 
state of the six- vertex model. This operator transforms under translations T x , T y as 

T x B{z)T x l = exp[-ip x (z)] B(z), (4) 
T v B(z)T- x = cxp[-z Py (z)]B(z), 

where equations 

exp[-ip x (z)] = Vk sn(z - (5) 



exp[— ipy(z)] = exp[— = Vk sn^z ■ , 

define in the parametric form the dispersion relation u>(p x ) of the one-particle 
excitations in the six- vertex model. Here sn is the Jacobian sn function to the modulus 
k, < k < 1, the quarter periods K, K 1 satisfy K'/K — 2X/tt. Parameters A and u 
are related with the original parameters of the Potts model by 0] : 

/ Q = 2 cosh A, (6) 
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exp(Ki) - 1 _ sinh[± A — it] 
~ sinh[± A + u] ' 
The elliptic modulus k can be explicitly written as 

1/2 

*- '-n(f 

n=l 

where 



• 1 — g 2 ™- 1 \ 8 



g2n-l 



(7) 




The one-particle spectrum (^|) can be obtained in the standard manner j?], |^| from 
the Bethe ansatz solution of the six- vertex model || . In the more general case of the 
eight-vertex model, such calculations have been done by Johnson, Krinsky and McCoy 
[0. It should be noted, that equations (||) describe also the one-particle excitation 
spectrum in the IRF and the critical Potts model (lj q|. In the critical Potts model 
these particles can be interpreted as the domain walls between an ordered and the 
disordered phases || . 

Let us transform the dispersion relation (|J) as follows. Consider transformation 
properties of the operator B(z) with respect to the translations along the rows and 
columns of the original Potts lattice. From ([|), (^) one obtains immediately 

T\ B{z)T{ 1 = expHPiW] B{z), (8) 
T 2 B{z)T^ = exp[-ip2(z)]B(z), 

where 

Pl(z) =Px{z) +Py(z), (9) 
Pz(z) = - p x (z) +p v (z). 

By use of (||), and applying the well-known properties of the Jacobian elliptic functions, 
one can exclude parameter z from (^|): 

ai cospi + 0,2 cosp2 = 1, (10) 

where 

ai = [(i + P 2 )(i + fcV)]- 1/2 , (ii) 

a 2 = fcp 2 [(i+ P 2 )(i + fcy)r 1/2 , (12) 

' 2 ^ + ^V (13) 
7T 2 / 

Since pi and e = ip 2 define the particle quasi-momentum and energy, respectively, 
equations ([l0|)-([l2|) give a compact form of the dispersion relation for one-particle 
excitations in the critical Potts model. The latter is exactly the same, as the excitation 
spectrum in the Ising model on the square lattice first obtained by Onsager jlO) . The 
only difference is in the meaning of parameters k and p. In the Potts model they are 
given by equations (Q), (|l3]), whereas in the Ising model we had 

k = [811111(2X1) sinh(2if 2 )] _1 , p = sinh(2if 2 ). (14) 

Substitution pi — > 1x2, p 2 — * ix\ into the dispersion relation ( [l0| ) yields |p"l| , |l2]| 
the equilibrium crystal shape (ECS) 

a\ cosh x 2 + 0,2 cosh x± = 1, (15) 
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where x%, x 2 denote the Cartesian coordinates in the equilibrium crystal boundary. 
This shape is the same in the critical Potts model [Q, and in the Ising model Q], if 
the temperature T < T c and the anisotropy ratio in the Ising model are determined 

by ©■ 

It is interesting to note, that the dispersion relation ( JlCj ) is relevant also to the 
simple Gaussian model on the square lattice, defined by the Hamiltonian 

P Eg = ^XX a( ^M + c i(f 3,1+1 - <Pj,lf + c 2 ((p J+ i,i - <Pj,i) 2 }, (16) 

where ipj 1 is a real continuous order parameter. Really, the two-point correlation 
function in this model can be written as 



where 



1 [ n dpi dp 2 exp(ipi I + ip 2 j) , . 

" 1 2 ci + 2 C2 J-tt {Sir) 1 - «i cospi - a 2 cosp 2 



2ci 2 c 2 

a l = > a 2 



2 ci + 2 c 2 a + 2 ci + 2 c 2 

The dispersion relation is determined by the pole of the integrand in (|l7|), i.e. by 
equation (|l0|). Since the one-particle dispersion relation determines the angular 
dependence of the correlation length, the latter is the same in a wide class of square- 
lattice models including the Ising, six-vertex, critical Potts, and the simple Gaussian 
model (|l6|). It is likely, that the origin of this universality lies in the symmetry of the 
square lattice. If this is the case, one could extrapolate this result to higher dimensions, 
speculating, for example, that the angular dependence of the correlation length in the 
three-dimensional Ising model on the cubic lattice could be identical with that in the 



Gaussian model like (16) on the same cubic lattice. 
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